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Abstract. We investigate the determination of the pole associated to a from Kit scattering data 
below the KK threshold, including the new precise data from K e 4 decay reported recently by the 
NA48/2 Collaboration. Using a large class of analytic parametrizations based on expansions in pow- 
ers of conformal variables, we obtain for the mass and width of a values which are consistent with 
those calculated recently using ChPT and Roy equations, but have larger theoretical uncertainties. 
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INTRODUCTION 

A precise determination of the mass and width of the o resonance using Chiral Pertur- 
bation Theory (ChPT) and Roy equations was obtained recently in [1]: 

M c7 = 441+ 1 8 6 MeV, T a /2 = 272+ 1 9 25 MeV. (1) 

In the derivation of this result, the isoscalar 5-wave is calculated at low energies, and 
also in the complex plane, from Roy equations, using experimental data at high energies 
and theoretical results on %% scattering [2, 3]. So, unlike in the standard way of detecting 
resonances, no experimental data on the partial wave with the quantum numbers of the 
resonance at low energies were used as input. The approach based on Roy equations is 
very suitable in this case, where the pole associated to the resonance is located far from 
the physical region and the experimental data at low energy are quite poor. 

Recently [4], NA48/2 Collaboration measured the phase shift difference 8$ — 5/ at 
low energies from K e 4 decay, with a precision much greater than that of the older 
experiments [5,6]. This revived the interest in the determination of the pole associated to 
<7 in the standard way, i.e. by the direct analytic extrapolation of the %% scattering data. 
In [7] the authors consider a representation of the isoscalar 5-wave t^{s) as an expansion 
in powers of a conformal mapping variable, and predict the mass and width of o with 
an accuracy comparable to that quoted in (1). 

In the present work we focus on the problem of systematic uncertainties within 
the approach proposed in [7]. By enlarging the class of admissible analytic functions 
considered in that work, we determine the mass and width of a from experimental data 
below the KK threshold, with a more realistic estimate of the uncertainties. 

Below we give a summary of the results obtained recently in this study. More details 
on the analysis and other results will be presented in a forthcoming publication [8]. 



ANALYTIC PARAMETRIZATIONS OF THE AMPLITUDE 



1. Writing the isoscalar S-wave t®(s) as: 



1 



P ( S ) = yJl-4M*/s, 



(2) 



elastic unitarity implies that is analytic in the s -plane cut only for s < and s > \M\ 
(we neglect the inelasticity due to the An channel below 1 GeV), except for a pole at 
the Adler zero sa, where ^(sa) = 0. The effective range approximation amounts to an 
expansion of \ir(s) in powers of s near s = AM 2 .. The domain of convergence can be 
enlarged by expanding in powers of a variable which conformally maps the holomorphy 
domain onto the interior of a disk [9]. The function 



w(s,a) = 



yfs-a^AM 2 -s 
y/s + aJ4M 2 -s 



(3) 



with a > arbitrary, maps the s-plane cut along s < and s > 4M| onto the unit disk 
\w\ < 1 in the complex plane w = w(s, a), such that w(4M|-, a) = 1 and w(0, a) = — 1. 
In [7] the authors adopt the expansion 



Y(s) = 



Mi 



SA 



_2sa_ 



+ B + B lW (s, a) + B 2 w(s, a) 2 + ... 



(4) 



with the particular choice a = 1 . In Eq.(4), the first term in parantheses compensates 
the singularity of p(s) at s = in the denominator of (2), removing a ghost of t®(s) on 
the real axis which would appear otherwise. A slightly different form was also used in 
[7], with an additional factor (/1q —s)/Hq inserted in (4). This parametrization, which 
displays the energy where the phase shift passes through 7t/2, is useful for fitting narrow 
resonances, but is not suitable for broad resonances like o, and we shall not use it. 

A first generalization of [7] is to expand in powers of w = w(s, a) with an 
arbitrary a, as in (4). By varying a, one changes the point mapped on the origin of the 
w-plane and the position of the intervals where experimental data are available. Some 
examples are shown in Fig. 1. 

2. An amplitude free of unphysical singularities can be alternatively obtained if the term 
ip (s) in (2) is replaced by a function which is analytic in the s-plane cut along s > AM\ 
and has the imaginary part equal to p (s) on the upper edge of the cut. We consider for 
convenience the loop function of ChPT: 



7(s,M2)) = I 



2 + p(s)ln 



Pis) -I 



(5) 



which vanishes at the origin, 7(0, M 2 n ) = 0, and has \mJ(s + i£,Ml) = p(s) for s > AM 2 % . 
Writing the partial wave as 



t° (s) = 



¥l ( s )-j(s,Miy 



(6) 



Re w 



Re w 



Re w 



FIGURE 1. The disk \w\ < 1 in the complex plane w = w(s, a) defined in (3), for a = 0.36 (left), a = 1 
(center) and a = 4 (right). The thick segments indicate the regions where experimental data are available 
from K e 4 decay [4]-[6] and the process 7tN — > %%N (cf. the compilation given in [10]), respectively; the 
circle shows the a pole on the second Riemann sheet from [1]. 



elastic unitarity implies that the function \\f\ (s) can be expanded most generally as 

= — — [Bo +Bi w(s, a) +B 2 w{s, a) 2 + ...]. (7) 

3. Other parametrizations of t^(s) are obtained by expressing the S-matrix element as a 
product 

S (s)=S rest (s)S fo (s), (8) 

where each factor satisfies elastic unitarity (|S r est(s)| = l^/o^)! = 1) below the KK 
threshold. It is convenient to parametrize the amplitude tf (s) associated to Sf (s) as 

tt ( s ) = ^ (9) 

/oW K-s-k lS J(s,Ml)-(k 2 + k 3 s)J(s,Miy 

where J(s, M 2 ) is defined in (5) and J(s,Mg) is obtained by replacing M n with Mk. We 
note that with the choice 

jc=1.01, fci=0.08, k 2 = -l.09, fc 3 = 1.16, (10) 

the modulus of S/ (s) above the KK threshold is close to the elasticity T]q(s) measured 
in [11], while for 

k= 1.15(1.41), ki = 0.11(0.24), k 2 = 0.39 (-0.73), k 3 = 0.03 (1.72), (11) 

it follows the upper (lower) edge of the experimental band of t/q (s) measured in [12]. 

In our fits, both the phase shift and the elasticity above the KK threshold are left free. 
The fact that we take a specific form for tf (s) is not a limitation, since t$(s) contains an 
additional term, t rest (s), related to the factor S rest (s) in (8). Elastic unitarity implies that 
^rest(s) can be parametrized, as in (2), in terms of a function y/r est (s) expanded as: 



/ N M l 

s — S\ 



— ^7= +B + B lW (s, a)+B 2 w(s, af + ... 



(12) 



where si is close to the Adler zero sa [8]. Admissible parametrizations are obtained also 
by replacing in this expansion w(s, a) by the variable w\ (s, a) = (y/s — a)/ (y/s + a), 
which maps the s-plane cut only for s < onto the unit disk \ w\ (*, a) | < 1 . 
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FIGURE 2. Left: phase shift 5q derived from K e 4 decay, fitted with the 16 parametrizations described 
in the text. Right: extrapolation of the parametrizations above the experimental range. 



RESULTS AND DISCUSSION 

We used the parametrizations described above for fitting the low energy data on the phase 
shift 8{f(s). In our analysis, the positive number a defining the conformal mappings 
and the parameters K and kt appearing in (9)-(ll) represent the input which defines an 
admissible class. In each class, the coefficients 5; of the expansion in powers of the 
conformal variable are free. They are determined by fitting the low energy data. 

We consider first the data on the difference 8q — 5/ measured below 0.4 GeV from K e 4 
decay [4]-[6]. The P-wave is known with precision in this energy range [3, 10], allowing 
an accurate extraction of phase shift of the 5-wave. As in [l]-[3], we work in the limit 
of exact isospin symmetry and take for convenience, for both M K and M K , the masses of 
the charged mesons. In order to obtain the strong phase shift 8q, an isospin correction 
calculated recently in [13] was subtracted from the measured phase shift. Following [7], 
we increased the experimental error on the last point in [6] by 50%. For the 10 data from 
the NA48/2 experiment we used the covariance matrix published recently in [4]. 

We investigated 16 admissible parametrizations: the first 3 are based on Eqs. (2) and 
(4) with a = 1, a = 0.36 and a = 4, respectively, the next 3 are based on Eqs. (6) and 
(7) with the same choices of a, and the last 10 parametrizations are based on Eqs. (8)- 
(12) with various conformal mappings [8]. In each case the Adler zero sa was allowed to 
vary between 0.4 M\ and 0.6 M\. Using 2 free parameters, Bq and B\, for 21 points, the 
optimal values of % 2 for the 16 parametrizations are: 21.7, 21.5, 21.9, 20.9, 21.2, 20.6, 
21.5, 21.7, 21.8, 21.6, 21.8, 21.5, 21.7, 21.8, 21.6, 21.4 (they decrease by about 0.4 units 
if the theoretical uncertainty of the isospin correction [13] is taken into account). 

The quality of the fits is seen in Fig. 2, where the experimental points are obtained 
from the data on K e 4 decay [4]-[6] as discussed above. Although the fits are almost indis- 
tinguishable in the experimental range, they exhibit large differences when extrapolated 
to higher energies, as seen in the right panel of Fig. 2. This illustrates the well-known 
instability of analytic extrapolation [14]. We note that the parametrization (6), especially 
with the choices a = 0.36 and a = 1 in (7), leads to phase shifts which exhibit a plateau 
at low values. The increase of <5q required by the high energy data is obtained, for in- 
stance, with a = 4 in Eq. (4), or by using the S-matrix factorization (8). 
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FIGURE 3. Positions of the a-pole obtained by the analytic extrapolation of the parametrizations fitting 
the various sets of data, compared with Refs. [1], [7] and [15] (from the last reference we show the value 
obtained with the isospin correction [13] included in the K e 4 data). 



By extrapolating to the second Riemann sheet of the s-plane, we obtain the positions 
of the a pole shown in Fig. 3, where m a = y/s^ = M a — iT a /2. The three isolated 
points, with small values for both the mass and width, correspond to the fits based on the 
parametrization (6) mentioned above. Taking the average of the 16 admissible values of 
M a and T a we obtain 

M a = 447±6(stat)+^(syst)MeV, r a /2 = 258 ±6(stat) +^(syst) MeV, (13) 

where the systematic error covers the spread of the admissible fits, including also an 
error of about 4 MeV in M a and 3 MeV in r,j, produced by the uncertainty in sa- 

We can improve the description of t® (s) by including data on the phase shift at higher 
energies. Two sets of data from kN — > kkN were considered: CERN-Munich data [16] 
consisting of 19 points below AM\, and a collection of 1 1 data points given in Eq. (2.13) 
of [10]. Using 13 parametrizations similar to those described above [8], we obtained, 
with 2 or 3 parameters 5,, values for x 2 m the range (33, 38) for the 40 points of the 
set I, and in the range (23, 29) for the 32 points of the set EL The quality of the fits is 
shown in Fig. 4, and the pole positions are given in Fig. 3. The three isolated points 
present in the fits of K e 4 data are no longer allowed. On the other hand, the narrow range 
of the widths T a , exhibited by the other fits of K e 4 data, is now enlarged. This is due to 
the fact that the description of the K e \ data, measured by their contribution to the total 
X 2 , is slightly worse than in the previous fits, and the various parametrizations are not as 
indistinguishable at low energies as in Fig. 2. 

Taking the average over the admissible parametrizations we obtain, for the two sets: 

M f7 = 455±6(stat)^3(syst)MeV, T a /2 = 277±6(stat)+^(syst)MeV (I) 

M a = 463±6(stat)+5 7 (syst)MeV, r a /2 = 259±6(stat)lg(syst)MeV (II). (14) 

Alternatively, we can define the central values by selecting the fits with the lowest values 
of x 2 hi each set. This procedure gives [8]: 

M (7 = 446±6(stat)+J°(syst)MeV, T a /2 = 267±6(stat)l^(syst)MeV (I) 
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FIGURE 4. Left: fits of data in set I (K e4 [4]-[6] plus CERN-Munich data [16] below the KK threshold). 
Right: fits of data in set II (K e4 plus a selection of data from nN — » MiN, given in [10]). 



Ma = 458±6(stat)+^(syst)MeV, T a /2 = 252±6(stat)l^(syst)MeV (II). (15) 

The comparison of (14) and (15) with (1) shows that the a-pole found by the analytic 
extrapolation of the low energy data on the kk isoscalar S-wave is consistent with the 
predictions of ChPT and Roy equations. However, the theoretical uncertainties are now 
larger, since the differences between the various parametrizations of the partial wave are 
amplified by the extrapolation from the physical region to a distant point in the complex 
plane. In the method based on Roy equations [1], the instability of the extrapolation is 
tamed by using additional information on the Tin amplitude. 
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